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Abstract 

A kinetic equation approach is applied to model anomalous J / ip suppression 
at RHIC and SPS by absorption in a hadron resonance gas which successfully 
describes statistical hadron production in both experiments. The puzzling 
Q-i| rapidity dependence of the PHENIX data is reproduced as a geometric effect 

r-| \ due to a longer absorption path for J/ip production at forward rapidity. 
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VO 

1. Introduction 

The effect of J/vp suppression [l[ is one of the diagnostic tools to prove 
G\ | quark-gluon plasma (QGP) formation in heavy-ion collisions. After a rich 

phenomenology provided by the CERN-SPS experiments NA3, NA38, NA50 
and NA60, which nevertheless could not unambiguously prove QGP forma- 
tion, the situation became even more complex due to first results from RHIC 
experiments PHENIX and STAR. New experiments are planned at CERN- 
LHC and FAIR-CBM, for a recent review, see j^j. In order to make progress 
in the interpretation of experimental findings it is important to define bench- 
marks, such as a baseline from a statistical kinetic model provided in this 
communication, which is based on a purely hadronic description. 

Measurements of J/ip suppression in Au-Au and Cu-Cu collisions by 
STAR and PHENIX experiments @, |, fl brought up two unexpected re- 
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suits. First, at forward rapidity J /if) is more suppressed than at midrapidity. 
Second, the dependence of J/ip suppression at midrapidity on the number 
of participants iVp art [1, 0] coincides with the one observed at SPS Pb-Pb 
collisions by the NA50 experiment [7|. Two different mechanisms to explain 
lata are usually considered in the literature: recombination 
or statistical coalescence 12, 13, 14, 15 and nuclear effects 
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We present a systematic step towards a general description of J/ip absorp- 
tion in the framework of a statistical analysis. The formalism allows 
for a unified description of the data from the NA38 and NA50 experiments at 
CERN and from the PHENIX experiment at RHIC. Both PHENIX results, 
concerning rapidity and N paTt dependences of J / ip suppression in Au-Au and 
Cu-Cu collisions are simultaneously reproduced. 

The J / ip absorption is caused here by the effective hadronic medium con- 
sisting of a mult i- component non- interacting hadron resonance gas (HRG). 
All hadrons from the lowest up to 2 GeV mass are taken into account as 
constituents of the matter |2l[. The gas is in thermal equilibrium and ex- 
pands longitudinally and transversally according to relativistic hydrodynam- 
ics [221 ] . J/ip suppression is the result of inelastic scattering on constituents 
of the HRG and on nucleons in the cold nuclear matter (CNM) of the collid- 
ing ions. Both sources of J/ip suppression, namely absorption in CNM and 
HRG, are considered simultaneously. 

The model applied here is a straightforward generalization of the model 
of Refs. 23j, ll9( to include the non-zero rapidity case. In our model, J/ip 
suppression is the result of a J/ip final state absorption in a confined medium 
through interactions of the type 



J/ip + h 



D + D + X, 



(1) 



where h denotes a hadron, D is a charm meson and X stands for a particle 
which assure conservation laws (charge, baryon number, strangeness). Ac- 
cording to the above assumption, charmonium states can be absorbed first 
in nuclear matter and soon after in the hadron gas. In our phenomenological 
analysis we assume universal cross sections for baryons, with the appropriate 
thresholds for their dissociation reactions but energy- independent, uj/^n = 4 
mb. This choice corresponds to the absorption cross section cr a b s = 4.2 ± 0.5 
mb for J/ip in cold nuclear matter measured in p-A collisions by the NA50 
experiment |24J , shown to be compatible 18|, [25] with recent d-Au data from 
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the PHENIX experiment |26|. The absorption cross section for meson impact 
is taken as a m = 2<tj/^ n /3, according to quark counting rules. 

An A-A collision at impact parameter b generates an almond shaped 
overlap region S e ff as presented in Fig. [TJ The time t — corresponds to the 
moment of the maximal overlap of the nuclei. After about half of the time the 
nuclei need to cross each other, matter appears in the central region. Soon 
thereafter matter thermalizes and this moment, t , is estimated at about 1 
fm. In the real situation the longitudinal thickness of the matter at to is 
also of the order of 1 fm. Then matter starts to expand and cool according 
to relativistic hydrodynamics (for details see 22, 19, [io[) until reaching the 
freeze-out temperature T f and subsequently freely streaming towards the 
detectors. 

A cc pair is created in a hard nucleon-nucleon collision during the passage 
of the colliding nuclei through each other. It will evolve to a charmonium 
eigenstate during the formation time which is of the order of to while experi- 
encing cold nuclear matter (CNM) effects. These are due to nuclear modifica- 
tion of the initial state gluon distributions and partial extinction while pass- 
ing through the colliding nuclei and will be appropriately parametrized using 
data from pA collisions. We are interested here to describe the "anomalous 
suppression" effect by subsequent final state interactions with the expanding 
HRG formed between the receding nuclei. Since these nuclei which border 
the HRG on both sides in the longitudinal direction move almost with the 
speed of light, even J/i/js which contribute to the particle spectra measured in 
the forward (backward) rapidity region can not escape the hadronic medium. 
Due to the different production process, the J / ip velocities can be considered 
as independent from the velocity distribution of the HRG medium. This re- 
sults in considerable impact velocities, in particular for forward (backward) 
produced J/ips traversing counter-propagating HRG flow. 

The transverse expansion starts in the form of a rarefaction wave moving 
inward S e ff. The evolution is decomposed into a longitudinal expansion in- 
side a slice bordered by the front of the rarefaction wave and the transverse 
expansion. The rarefaction wave moves radially inward with the sound ve- 
locity c s . Since the temperature (and hadron gas density) decreases rapidly 
outside the wave, we shall ignore possible J/ip scattering there. We denote 
by t esc the moment of crossing the border of the rarefaction front. 

Inside the region bordered by the front of the rarefaction wave the hydro- 
dynamic evolution ceases when the freeze-out temperature is reached, here 
we take Tf. Q = 150 MeV what is suggested by the statistical model analy- 
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Figure 1: View of a A- A collision at impact parameter b in the transverse plane s = 
(x, y, z = 0). The region where the nuclei overlap has been hatched and its area is S e ff. 



sis of the PHENIX data [27|, |28|]. So J/ijj is subject to absorption for the 
time tfinai = min{(t esc ), tf.o.} where (t esc ) is the mean time until it escapes 
the hadron gas region and tf. Q . is the time when it passes the kinetic freeze- 
out surface for its absorption reactions. The rapidity dependence of tg na i is 
crucial for the description of J/if> RHIC absorption processes. 

For comparison with data one should estimate the number of participating 
nucleons N part (b). The Woods-Saxon nuclear matter density distribution 



for Au and Cu nuclei is assumed with parameters taken from [29| and the 



standard expression for A^ part (6) is used |30l. Il9 



2. J /if) absorption 

The J /if) absorption processes in CNM (nuclear absorption - NA) and 
HRG take place subsequently. Due to this separation in time the J/if> survival 
probability for given rapidity in a heavy-ion collision with impact parameter 
b assumes the form 

Raa(v, b ) = S NA (y, b) ■ S nRG (y, b) , (2) 
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Figure 2: J/i/> (dashed and dot-dashed lines) inside hadronic gas absorption region and 
the time evolution of the rarefaction front radius s at different z (continuous lines with 
labels denoting z in fm). Numerical values correspond to Au-Au collisions. 



where Sna(2/> b) and <Shrg(z/>&) are J/ip survival probabilities in CNM and 
HRG, respectively. For Sna we employ the usual approximation 

S'na = exp{-a^ N p (L)} , (3) 

where (L) is the mean absorption path length of the J/ip through the colliding 
nuclei obtained from the Glauber model and po is the average nuclear matter 
density determined from the normalization of the Woods-Saxon distribution 
to the mass number A for each nucleus separately. Within this approximate 
expression, S'na does not depend on rapidity. 

The HRG survival probability Shrg(z/> b) is defined as 

/ d 2 p T J d 3 r F(r,y,p T ,t)\ 
J d 2 p T J d 3 r F(r, y, p T , t) \ t =t 

where the initial time to = 1 f m denotes the moment of the thermalization 
of the created hadronic matter and the beginning of the hydrodynamical 
expansion. 

The rapidity- momentum J/ip distribution results from the kinematical 
change of variables d 3 p —> d 2 pxdy and is given by J-(r, y,pr, t) = M T cosh(j/) 
xf(f,p, t), where 



HRG (l/j b) = -jr— ^— , (4) 



f(r,p,t) = fo(r-v(t-t ),p) 



(5) 



x exp 



/'*'£ / 

Jto i=1 J 



d 3 q i PuQ: 

fi(r -v{t-t),q,t)aiV rd 



(27T) 



EE' 



is the formal solution of a kinetic equation (for details see 23|, |l9|). The sum 



is over all considered species of hadronic scatters with distributions /j(r, q, t). 

We have assumed that the initial distribution factorizes into J-(r, y, px, to) 
— foif) ' 9o(Pt) " ho(y) with fo(f) and <7o(pt) normalized to unity. It is 
clear that /i (j/) has not to be specified since it cancels in the ratio (jl]). 
Because of the arguments given in the Introduction we shall assume that 
both the created J/ip and the hadronic medium are included in z = plane. 
Thus the distribution of J/ip at t is fo(s) where s belongs to the z = 
plane. For qniv r) we take the form suggested by the initial state interactions 



2p T 



Pt 



(p 2 t)^ 



(6) 



where (py)jA is the mean squared transverse momentum of J/ip measured 
in an Au-Au (Cu-Cu) collision 

For the description of the evolution of the matter, relativistic hydrody- 
namics is employed. The longitudinal component of the solution of the hy- 
dro dynamic equations (the exact analytic solution for an (l+l)-dimensional 
case) reads |35| 

s o r o 



ST 



n B {T 



n° B r 

T 



(7) 



where r = y/j? 



is the proper time and so (n B ) is the initial density of 



entropy (baryon number). 

For ub = and a uniform initial temperature distribution with a sharp 
edge at the border established by the nuclear surfaces, the full solution of 



the (3+l)-dimensional hydrodynamic equations is known 22 
The hadron distributions have the usual form: 



fi(r,q,t) 



2 Sl + l 



cxp 



q»Uy(r,t)-Hi(r,t) 
T(r,t) 



(8) 



+ 9i 
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where u v (r,t) is the four-vector the hadron flow velocity and B iy Si, 
Hi, Si and gi are the mass, baryon number, strangeness, chemical poten- 
tial, spin and a statistical factor of species i, respectively. We treat an an- 
tiparticle as a different species and introduce the chemical potential profile 
fj>i(f,t) = Bip,B(r,t) + Sifis(r,t), where /xs(r, t) and /xs(r,t) are the local 
baryon number and strangeness chemical potentials, respectively. 

Since for the flow described in [22] the radial velocity is zero inside the 
region bordered by the rarefaction wave and v z (s, z, t) = z/t, the four- velocity 
u v {f,t) simplifies to 

u»(r,t) = -(t,0,0,z) . (9) 
T 

Note that a J/ip with longitudinal velocity Vl is at z = vi(t — t ) at the 
moment t and the longitudinal flow at this point equals v z (s,z,t) = z/t = 
VL(t—to)/t < vl, that is the relative longitudinal velocity of J/ip with respect 
to the flow is always non-zero. This reflects the fact that only the evolution 
of the medium is boost-invariant but not of the J/ip. The latter is moving 
along a straight line with the constant velocity from the point of creation as 
it does not take part in the thermalization of the hadronic medium. 

The temperature also does not depend on the radial coordinate s within 
the region bordered by the rarefaction wave, 

T(f, t) = T{z, t) = T(0, r) = T ■ r~ a (10) 

and the power a is the sound velocity squared at T , a = c^(T ) 0, 20]. 



To obtain To and a the following procedure is applied. For a given cen 
trality bin the initial energy density eo is estimated in (Au-Au) and 



(Cu-Cu). We put ns = for the strangeness number density and n° B = 
for the initial baryon number density. It should be noted here that the J/ip 
suppression pattern is practically the same for values of n° B < 0.65 fm -3 
in this model as it was shown for the case of the centrality dependence in 



Ref. [19j]. This range of n° B corresponds to the possible values of [lb up to 



300 MeV during the whole evolution. This is far above the estimated val- 



ues of hb = 20 — 30 MeV at the top RHIC energy [28|, [38|. Now, expressing 
e , n° B and ns as functions of T, \i B and \i$ in the grand canonical ensemble 
framework, we can obtain T , n B and (j? s . In our case fiB = H°b = an d 
fis = A*s = always, so the temperature is the only significant statistical 
parameter here. 

For the most central collisions we have obtained T = 222.6 MeV (eo = 5.4 
GeV/fm 3 ) for Au-Au and T = 201.8 MeV (e = 2.5 GeV/fm 3 ) for Cu-Cu. 
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With T known the initial entropy density s is calculated. Having put s into 
(J7j) and with the use of the grand canonical expressions for s we can obtain 
T(t). It turned out that T(r) has the form given by (fit)]) with a varying in 
the range a = 0.148 — 0.156 for eo in the range eo = 5.5 — 0.5 GeV/fm 3 . 

For a J/ij} which is at sq G S e ff at the moment to and has the velocity 
v = Vt +Vl = Pt/E + Pl/E we denote by t esc the moment of crossing the 
border of the rarefaction front. 



The rarefaction front forms the ellipse in the z, s plane [22 

(s - R A ) 2 + c 2 s z 2 = c 2 s t 2 (11) 

The continuous lines in Fig. [2j labeled by different values of the longitudinal 
variable z, show the time evolution of the rarefaction front radius s at the 
position z. 

The J /if) trajectory is given by equations 

8 = d + v T (t-t ), (12) 
z = VL (t-to). (13) 

We have used here the notation d = s — b for the J/ ifi which escapes via the 
left half of S e ff and d = s in the opposite case. 

The dashed and dash-dotted lines in Fig. [2] are (s,t) trajectories of J /if) 
particles originating from the same point in the region S e ff shown in Fig. 1, 
but with different values of transverse velocity vt- The dot-dashed line cor- 
responds here to the lower value of vt- 

The rarefaction front crossing time t esc is a solution of the equation 



^d 2 + 2d- v T (t - t ) + v 2 T (t - t ) 2 = R A - c s ^t 2 - v 2 L (t - t ) 2 . (14) 

Having obtained t esc , we average it over the angle between s*o and vt- 
Finally, to obtain {t esc )(b,VT,VL) - the average time it takes the J /if) to leave 
the hadronic medium, when it is produced in a collision at impact parameter 
b with the velocity = s/v 2 + v 2 and Vl - we average the result over S e ff 
with the weight given by the initial J /if) distribution /o(so)- This weight is 
equal to 

T A (s )T B (so - b) 

m«o) = — — ( 15 ) 

J-abK.0) 
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where 




Ta(s) — Jdz Pa(s, z) is the nuclear thickness function and pa(s, z) the Woods- 
Saxon nuclear density distribution. 

Then the final expression for Shrg(2/, b) reads 



SHR G (y,b) = -j- — — - I dp T M T g (p T ) (16) 

x exp | - dtJ2J J^fc&^^i^j > 

where tfi na i = min{(t esc ), t/. .} and t/. D . is the freeze-out moment resulting 
from the longitudinal Bjorken expansion. 

Eq. ( IT6|) is in fact an approximation of the more involved exact formula. 
The approximation means that the average of integrals over trajectories has 
been replaced by the integration over the averaged trajectory, similarly as in 
Eq. 03]). Our preliminary results of estimations of the exact formula suggest 
that the approximation influences only the normalization and in such a way 
that this can be compensated by the rescaling of the absorption cross-section 
to the lower values. 



3. Results 

In Fig. [3] an example of the behavior of tfi na i with rapidity is presented. 
This is for the case of the — 20% centrality bin of Au-Au and Cu-Cu col- 
lisions and for px = 2 GeV, taken as a typical value for collisions with a 
(p2) w ^4GeV 2 0,3. 

In Figs. |U <^ we show the J /if) nuclear modification factor Raa for 
different centralities versus rapidity in Au-Au (Cu-Cu) collisions vs. PHENIX 
data. There is an overall agreement with the data. 

The main new effect of the PHENIX data, that J /if} suppression is stronger 
at forward rapidity than at midrapidity is shown in Fig. [6] for the iVp ar t de- 
pendence. It arises in the present approach as a result of the geometry of 
the kinetic freeze-out surface for the J/if> which extends much more in the 
longitudinal direction for sufficiently long times, t > Ra — 6/2. This entails 
a longer absorption path for J /if) production at forward rapidity. 



9 




Figure 3: An example of ta na i in a function of rapidity for px — 2 GeV and the — 20% 
centrality bin for Au-Au (solid) and Cu-Cu (dashed) collisions. 



Note that the present model simultaneously accounts for the anomalous 
J/i/j suppression in the centrality dependence of the NA50/NA60 experi- 



ments at CERN, see [19l . |20j . Can this result, obtained within a purely 
hadronic description, be interpreted as evidence against QGP formation in 
these experiments? Not at all, because of two reasons. 

First, we have employed hadronic absorption cross sections which cannot 
be reconciled with a microscopic description in, e.g 
or relativistic 



39. 40 



41, 42 



nonrelativistic 
quark models which have a fastly decreasing energy 

In order to 



dependence and do not exceed a peak value of 2 mb, see also [43 



justify the magnitude of hadronic absorption cross sections employed in the 
present work, a strong medium dependence is required which could for in- 
stance stem from spectral broade ning of light mesons at the chiral restoration 
transition (mesonic Mott effect) |44j . 

Second, the hadronic resonance gas model gives a perfect decription of the 
Lattice QCD thermodynamics in the vicinity of the critical temperature T c 



for the chiral and deconfmement transition 45, 20 . It has been demonstrated 



that the description of QCD thermodynamics within a hadronic basis can be 
extended even up to temperatures of 1.5 T c , provided a suitable spectral 
broadening of these states above T c due to their Mott effect is taken into 
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Figure 4: J / ip nuclear modification factor versus rapidity in Au- Au collisions for constant 
Gj/ipN — 4 nib. PHENIX data are from Q. Errors shown are the quadratic sum of 
statistical and uncorrelated systematic uncertainties. 



account [46(. For a recent development, see |47|. Such a Mott-Hagedorn 
resonance gas model has been used to describe anomalous J /ip suppression 
at SPS jigj and is in accordance with the picture of a strongly coupled QGP 
(sQGP) discovered at RHIC 0, M ■ 

Summarizing this discussion, the present model may be seen as a prepara- 
tory step towards a unified description of charmonium suppression kinetics 
in SPS and RHIC experiments within a quantum statistical approach to the 
sQGP to be developed. The geometrical effect on the rapidity dependence 
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Figure 5: Same as Fig. |4]but for Cu-Cu collisions with PHENIX data from 
decribed in the present work shall be a part of such a description. 
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